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Abstract 

The thermodynamic behavior of QCD matter at high temperature is cur- 
rently studied by lattice QCD theory. The main features are the fast rise of 
the energy density e around the critical temperature T c and the large trace 
anomaly of the energy momentum tensor (0^) = e — 3P which hints at a 
strongly interacting system. Such features can be accounted for by employ- 
ing a massive quasi-particle model with a temperature-dependent bag con- 
stant. Recent lattice QCD calculations with physical quark masses by the 
Wuppertal-Budapest group show a slower increase of e and a smaller (0^) 
peak with respect to previous results from the hotQCD collaboration. We 
investigate the implications of such differences from the point of view of a 
quasi-particle model, also discussing light and strange quark number suscep- 
tibilities. Furthermore, we predict the impact of these discrepancies on the 
temperature-dependence of the transport properties of matter, like the shear 
and bulk viscosities. 

PACS: 12.38.Mh, 25.75-q, 52.25.Fi 

1 Introduction 

QCD at very high temperature T is expected to be weakly coupled pQ and provides 
perturbative screening of the charge [2J, thus being called Quark-Gluon Plasma 
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(QGP). Creating and studying this phase of matter in the laboratory has been 
the goal of experiments at CERN SPS and at the Relativistic Heavy Ion Collider 
(RHIC) facility in Brookhaven National Laboratory, recently continued by the AL- 
ICE, ATLAS and CMS Collaborations at the Large Hadron Collider (LHC). These 
experiments provide the unique possibility of quantifying the properties of the decon- 
fined phase of QCD. At the same time, lattice calculations on QCD thermodynamics 
are now reaching unprecedented levels of accuracy, with simulations at the physical 
quark masses and several values of the lattice cutoff, which allows to keep lattice ar- 
tifacts under control. The information that can be obtained from these two sources 
can shed light on the features of QCD matter under extreme conditions, one of the 
major challenges of the physics of strong interactions. 

As soon as the first lattice results on the QCD equation of state became avail- 
able, several attempts have been made in order to interpret them in terms of the 
appropriate effective degrees of freedom. 

There have been many efforts to discuss QCD thermodynamics at high tem- 
peratures by means of pertubation theory. For the pressure, the calculations have 
been extended up to g 6 log g [3] , but besides strong fluctuations from one order to 
the other, the final result agrees with lattice data only for T > 5T C . Moreover, if 
one concentrates on the interaction effects by calculating the interaction measure 
(6p = e — 3P, the pertubative expansion completely fails to reproduce the lattice 
data. In order to overcome this problem, resummation schemes have been proposed, 
based for example on the Hard Thermal Loop (HTL) approach [3] |5| 16] IT] 18] 19] [TO] [TT] 
or on dimensionally reduced screened perturbation theory (DRSPT) [J31 HU [J3] . 

More recently within a NNLO in the HTL perturbation theory [16J it has been 
shown that an agrement with the lattice results can be obtained also for the trace 
anomaly down to T ~ 2T C . At these high temperatures, the HTL approach moti- 
vates and justifies a picture of weakly interacting quasi-particles, as determined by 
the HTL propagators. 

A particularly intriguing possibility is to bring the description of lattice results 
in terms of quark and gluon quasi-particles down to T c . Their effective masses are 
generated through the interaction among the basic constituents [T71 HH HH 120] : 
if a large part of the interaction can be included into the effective masses, quasi- 
particles move freely or interact only weakly. This approach can be very useful in 
phenomenological hadronization models [2JJ and also for a description of the QCD 
medium in the vicinity of the phase transition, where perturbative methods cannot 
be used. The quasi-particle model (QPM) was successfully applied to quantitatively 
describe lattice QCD results of equilibrium thermodynamics such as the EoS and 
related quantities [22]: by a fit to lattice QCD thermodynamics, the quasi-particle 
properties such as their effective masses can be obtained. 

Recently, new lattice results for the equation of state of QCD with 2+1 dynamical 
flavors have become available [23J. These simulations have been performed at the 
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physical quark masses and for several lattice spacings, and a continuum estimate of 
the bulk thermodynamic quantities has been given. These new results show several 
differences with respect to the previously available ones [2U [25] (a slower increase 
of e and a smaller (0^) peak), which stimulates the question whether and how the 
properties of quark and gluon quasi-particles are affected by these features. 

In the present paper we assume that it is possible to describe the QCD medium 
above the deconfmement phase transition in terms of quark and gluon quasi-particles 
with an effective mass. Our purpose is to check how the properties of these quasi- 
particles change, following the differences between the available lattice data for QCD 
thermodynamics. After performing a fit to the available lattice data, we will obtain 
the quark and gluon temperature-dependent masses, which are related to each other 
through the coupling constant, and a temperature-dependent bag constant which 
is necessary for thermodynamic consistency. Once the thermodynamic properties 
of the medium are determined, we will extract dynamical properties such as bulk 
and shear viscosities and relaxation time. The model predictions for quark-number 
susceptibilities are also discussed. 

2 The Quasi-Particle Model 

A successful way to account for non-perturbative dynamics is a quasi-particle ap- 
proach, in which the interaction is encoded in the quasi-particle masses. The mass 
of the particles can be viewed as arising from the energy contained in a strongly 
coupled volume determined by the correlation range of the interaction. Once the 
interaction is accounted for in this way, the quasi-particles behave like a free gas of 
massive constituents. The model is usually completed by introducing a finite bag 
pressure that can account for further non-perturbative effects and could be directly 
linked to the gluon condensate at least in the pure gauge case It is already 

very well known that, in order to be able to describe the main features of lattice 
QCD thermodynamics, a temperature-dependent mass has to be considered. This 
also implies that the bag constant has to be temperature-dependent, in order to 
ensure thermodynamic consistency. In fact, when a temperature-dependent mass is 
included in the pressure, its derivative with respect to the temperature will produce 
an extra term in the energy density, which does not have the ideal gas form. A 
solution to this problem has been proposed in [T7]: pressure and energy density 
contain additional medium contributions, which we will call B(T), and which cure 
the problem of thermodynamic consistency, as we will see in the following. In this 
formulation, the entropy density preserves the ideal gas form. 

The temperature-dependent effective mass for quarks and gluons can be evalu- 
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ated in a perturbative approach that suggests the following relations [TU [19j 120] : 



m. 
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where n/ is the number of flavors considered (2 or 3 in this work, see the following), 
N c is the number of colors, m u ^ is the mass of the light quarks and \x q is the chemical 
potential of the generic flavor q considered. 

The coupling g is generally temperature-dependent. However, as mentioned in 
the introduction, the calculation of such a T-dependence by means of perturbation 
theory does not allow to have a good description of lattice QCD thermodynamics. 
Therefore, usually g(T) is left as a function to be determined through the fit to 
lattice QCD data. This will be also our methodology, but we will consider lattice 
QCD data for SU (3) in flavor space, hence for the strange quark mass we will use 
the following relation: 
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where m s is the total mass of the strange quark while mo s is the current quark mass. 
We see that Eq. ^ reduces to the second equation in (jl]) when m 0s = 0. The effect 
of a finite current quark mass is generally negligible for the description of the energy 
density, pressure, entropy but can play a role in the evaluation of the strange quark 
number susceptibility, see Section |4j 

The pressure of the system can then be written as the sum of independent contri- 
butions coming from the different constituents, which have a T-dependent effective 
mass, plus a bag constant: 

(2tt)3 3E-(p) ^ ~ B(yT ^ ' ® 

i=u,a,s,g 

where fi(p) = [1 =F exp are the Bose and Fermi distribution functions, 

with Ei(p) = a/p 2 + mf ; di = 2 x 2 x N c for quarks and di — 2 x (N^ — 1) for 
gluons. 

In order to have thermodynamic consistency, the following relationship has to 
be satisfied: 

^) =0, i = u,d,... , (4) 
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which gives rise to a set of equations of the form 
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Only one of the above equations is independent, since the masses of the constituents 
all depend on the coupling g through relationships of the form: mj(T, /i = 0) = 
aig(T)T, where «j are constants depending on N c and Nf according to Eqs. ([I]). 
The energy density of the system is then obtained from the pressure through the 
thermodynamic relationship e(T) = TdP(T)/dT — P(T) and will have the form 
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m t ,T) + B( mi (T)). (6) 



In the model there are therefore two unknown functions, g(T) and B(T), but they 
are not independent, they are related through the thermodynamic consistency rela- 
tionship Q. Therefore, only one function needs to be determined, which we do by 
imposing the condition: 

e qp (T) = ei a ttice(T). (7) 

Following Ref. [29J, we derive both sides of the above equation with respect to T, 
then by using equation ^ and recalling that m,i(T) = ceig(T)T, we obtain: 
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where 6km = ^2i e kin- Therefore, by performing the fit to the lattice energy density, 
we can solve the differential equation for g(T), given the initial value g(T ) = g . 
The function g(T) obtained in this approach for T > T c is the same as the one 
obtained in other models, in which a parameterized form for g(T) is assumed, and 
the parameters are fixed through a fit to the lattice data (see for example |20j): 
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In our model, g(T) is obtained numerically from Eq.([8]). If we want to express it 
in terms of the above parameterization (91), we obtain A = 17.7 and Ts/Tc = 1.15 
from a fit to the results of the Ref. |25r and A = 2.6 and Ts/Tc = 0.57 from a fit 
to the results of the Wuppertal-Budapest collaboration |23j. 

1 Notice that preliminary results for some observables have been published by the hotQCD 
collaboration with the new hisq action [27, 28 , which has a small taste violation at low temperatures 
and small cutoff effects at large temperatures. However, for our analysis we use the results obtained 
with the p4 action, since the energy density data are not yet available for the hisq action. 



5 




Figure 1: Lattice data for the pressure and trace anomaly as functions of the temper- 
ature, together with the quasi-particle model curves. The results of the Wuppertal- 
Budapest collaboration are from Ref. (23], the other lattice results are taken from 
Ref. [2S] • 



We recall here that the values of T c are different in the two sets of lattice data, 
more specifically we have used T c = 175 MeV for the hotQCD case and 155 MeV 
for the WB one. 

We perform our fit to the lattice data for the energy density; in Fig. [T] we show 
the good agreement between our curves and lattice results for other quantities like 
pressure and trace anomaly, due to thermodynamic consistency. 

In Fig. [2]we show the temperature dependence of the coupling g and of the quark 
and gluon masses that we obtain from this procedure. As we can see, the slower 
increase of the energy density and the smaller values at hight T in the Wuppertal- 
Budapest (WB) data with respect to the hotQCD ones is reflected into larger quark 
and gluon masses. At small temperatures, as T — > T c , the energy density from 
the Wuppertal-Budapest collaboration is larger than the hotQCD one, and as a 
consequence in the last case quark and gluon masses have a stronger T dependence 
resulting in larger masses in this temperature regime. 

We notice that, at sufficiently high temperature, m ~ T, as we can expect 
because T remains the only scale of the problem. When we are approaching the 
phase transition, there is instead a tendency to increase the correlation length of the 
interaction: the quasiparticle model tells us that this can be described as a plasma 
of particles with larger masses. This essentially determines the fall and rise behavior 
of m(T) seen in all quasiparticle model fits to lattice data in SU(N) gauge theories 
including the SU(3) case for QCD. In Fig. [2]we can see a quite smoother behavior of 
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Figure 2: Left panel: temperature-dependent coupling g(T) as a function of T/Tc- 
Right panel: quark and gluon quasi-p article masses as functions of T/Tc- We used 
iriQ S = in the fit: the effect of a finite current quark mass for the strange quark is 
negligible for the description of thermodynamic quantities. In both panels, the red 
curves correspond to a fit of the lattice results from the hotQCD collaboration |25J, 
the black ones to a fit of the Wuppertal-Budapest collaboration results [23J . 



the masses associated to the Wuppertal-Budapest case, indicating that the strength 
of such correlation is significantly reduced when lattice simulations are performed 
at the physical quark masses and the continuum limit is taken. 

The function B(T) can be obtained from the difference between the lattice energy 
density and the kinetic contribution from our model: B{T) = €i a ttice(T) — e kin (T). 
We show its temperature dependence in Fig. [3] for the two different fits to the lattice 
results. B(T) provides another measure for non-perturbative physics which cannot 
be absorbed into effective quark and gluon masses. As it is evident from Fig. [3j in the 
vicinity of the phase transition the B(T) corresponding to the hotQCD lattice data 
(red curve) is much larger than the Wuppertal-Budapest one (black curve). This 
mainly reflects the discrepancy in the interaction measure (e — 3P)/T 4 between the 
results of the two collaborations (see Fig. 17 of Ref. [23]) . ^From the energy density 
and the pressure one can easily evaluate the speed of sound squared, defined as 

d = f ; (10) 

we show it for completeness in Fig. |4j The curve corresponding to the hotQCD fit 
displays a steeper rise in the vicinity of the transition temperature, and approaches 
the Stefan-Boltzmann limit much faster with respect to the curve corresponding to 
the Wuppertal-Budapest fit. For comparison, in Fig. [4] we show the speed of sound 
as calculated for the Nambu-Jona-Lasinio (NJL) model, the set of parameters used 
are those taken by Buballa [30, 32J. The speed of sound calculated in the NJL model 
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Figure 3: B(T) as a function of T/Tq- 



has a rapid increase near the critical temperature and for T « 1.27c it converges 
towards the characteristic value of a relativistic massless gas. This is due to the fact 
that in this model the trace anomaly is very small for T > Tc, see [32]. However 
if the Polyakov loop dynamics is added the behavior of the sound velocity c s (T) is 
very well reproduced as shown by the long dashed line taken from Ref. |31J . which 
is fitted to the hotQCD data. 

3 Transport coefficients 

Once the properties of quark and gluon quasi-particles are fixed to reproduce the 
lattice data, we can use the model to compute quantities that cannot be evaluated 
on the lattice, such as for example transport coefficients. To study their thermal 
properties we use the method described by C. Sasaki and K. Redlich [33], which is 
based on the relativistic kinetic theory formulated in the relaxation time approx- 
imation. The formulas for the viscosities r\ and ( are derived for a quasi-particle 
description with bosonic and fermionic constituents (see [33J for details): 
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Figure 4: Speed of sound squared as a function of T/Tc- The red curve corresponds 
to the hotQCD data, the black one to the Wuppertal-Budapest data, the short- 
dashed one to the NJL model and the long-dashed one to the PNJL model fit to the 
hotQCD lattice data. 



In the relaxation time approximation, both shear and bulk viscosities for the parti- 
cles i — u,d, s, g depend on their collisional relaxation time Tj. The latter is given by 
the thermal average of the in-medium cross section a (describing total elastic scat- 
tering of medium constituents) times the relative velocity of two colliding particles 

Vrel [31]: 

r t r\T)=n t (T)(v rel a)(T); (13) 

rii is the density of the particle i. The in- medium cross sections for quark-antiquark, 
quark-quark and antiquark-antiquark scattering processes were studied in detail in 
Ref. [35] within the NJL model for two different flavors, including 1/Nq next-to- 
leading order corrections. These results incorporate dominance of the scattering on 
large angles and take into account a possible occupation of particles in the final 
state, see Ref. [36J. 

The QCD calculations of the relaxation time r of partons require the sum of 
infinitely many diagrams already in the lowest order in the running coupling constant 

9- 

A first evaluation of the hard thermal loops resummation [371 EH] has suggested 
a parton width r _1 going as ~ g 2 T ln(l / g) . Thus in recent papers [19], [39J HO] within 
a quasiparticle approach the following parametrizations have been used for quarks 
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and gluons: 
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where g is the coupling obtained by solving Eq. (|8j) and is a parameter which is 
fixed by requiring that 7$ yields a minimum of one for the quantity Airrj/s. 

We point out that such expressions do not recover the leading-order pQCD limit 
where r is proportional to g^lnfy -1 ) [13] and not g 2 ln(g~ 1 ). However, a pure g A 
parametric dependence is strictly valid in the weak coupling limit for T >> Aqcd, 
a condition that is not fully satisfied in our temperature range. Nevertheless in 
Ref. [44J it has been shown that parametric corrections to the leading order are 
subdominant and give at most a g power correction at least when mu < T. One 
may still argue that in our temperature range, especially at T < 2T C one is far 
from a weak coupling limit and the last condition is not fulfilled. Nevertheless we 
will consider the relaxation time expressions in the weak coupling limit showing 
that in the asymptotic limit our expression for rj/s can recover the correct pQCD 
limit. Indeed, considering both the g 2 and g A dependences for r _1 will be particularly 
interesting when more accurate 1QCD calculations for the shear viscosity will become 
available, allowing to infer the parametric dependence of r within a quasiparticle 
approach. We have considered the following ansatz for the relaxation time as in 
Refs. [HH HI] derived in a pQCD approach: 



t; 1 = C g g*T\n{a 2 /g 2 ) 

r- 1 = C q g*T In (a 2 /g 2 ). (15) 

where C q /C g = (JVg - 1)/(2N$) = 4/9. The validity region of pQCD is certainly 
not clear, however we will show in the following that the g A pQCD parametric de- 
pendence of the relaxation time gives a significantly different behavior as a function 
of the temperature. We use these Tj in order to recover, at large T, the dependence 
of the shear viscosity on temperature and coupling known from pQCD |43| |4*4] . 

In Fig. [Bl we show the relaxation time for gluons as a function of T/Tc for the 



two different lattice fits and the two different relaxation times in Eq.(14) and Eq.(15) 
(discussed later in more detail). The results shown in the figure are obtained using 
the same value of the parameter, k = 22, similarly to (42]. As we can see for the 
g 2 parametric dependence (dashed lines) the relaxation time is a monotonically de- 
creasing function of T for both lattice fits, apart from a shallow minimum slightly 
above T c in the hotQCD case. The trend changes significantly for the g A dependence, 
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Figure 5: Relaxation time r g for gluons as a function of T/Tc- The red curves 
correspond to the hotQCD data, the black curves to the Wuppertal-Budapest ones. 



The solid lines correspond to Eq. (14), the dashed lines to Eq. (15). 



that in both cases generates a minimum in the relaxation time: this can be con- 
sidered as an indication of a strong increase of the interaction just above the phase 
transition temperature. We notice that the smoother phase transition dynamics of 
the WB lattice QCD data is reflected also in the relaxation time, by generating a 
less deep minimum with respect to the hotQCD case. 

In Fig. [6]we show the shear viscosity to entropy density ratio, r//s, as a function 
of T/Tc for the g 2 case (solid lines) and the g 4 case (dashed lines). It is interesting 
to notice that, once k is fixed to have a minimum rj/s = 1/47T, this minimum 
appears very close to T c and is followed by a divergency as T — > T c . This behavior is 
shared by all the cases considered. Of course the exact location of the minimum goes 
beyond the parton quasiparticle approach, that breaks its vality for T < T c , however 
on general arguments a minimum is expected very close to T c jl5] , furthermore one 
also expects a large increase of T]/s below T c , entering the hadronic phase. 

At high temperature the viscosity corresponding to the Wuppertal-Budapest fit 
shows a slower increase with temperature with respect to the one corresponding to 
the hotQCD fit. In the case of the g 2 dependence of r, the increase in r]/s is quite 
mild and it seems that even at temperatures T ~ 4T C this transport coefficient is 
much smaller than its estimate in pQCD, rj/s ~ 1. This is in apparent contradiction 
with the behavior of the trace anomaly, which decreases quite rapidly with increasing 
T, approaching the perturbative estimate relatively early. Such a behavior can be 
mostly related to the uncertainty in the evaluation of the relaxation time r and in 



particular its temperature dependence. In fact, as mentioned before, Eq.(14) does 



11 



0,1 



- 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


- 














WB, x'-g'TlnCg') 








hotQCD 








WB, x" 1 ~g 4 Tln(g 1 ) 








hotQCD 


















/ , ----- 












I i ' 


^ — 

^ ✓ 






\y\z> — 




1 


r|/s=(47i)" 1 - 
I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 



1 1,5 2 2,5 3 3,5 4 
T/T c 



Figure 6: Shear viscosity to entropy density ratio rj/s as a function of T/Tc- The red 
curves correspond to the hotQCD data, the black curves to the Wuppertal-Budapest 
ones. The solid lines correspond to a g 2 dependence of r, the dashed lines to a g 4 
dependence of r. The dotted line is the rj/s — l/(47r) bound. 



not give the proper pQCD limit at very large temperature and may be questionable 
already at moderate temperature T ~ 4T C as those considered here. By using 



Eq.(|15|) not only the pQCD relaxation time is recovered, but also the correct limit 

In 

3 and N c 



for the shear viscosity is guaranteed by Eq. (11), as we discuss in the following. 
jH] the shear viscosity was obtained for Nf 
order of the small running coupling g. 
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with rji = 106.66 and /x*/T = 2.957. In our quasiparticle approach, from Eq.(ll) 



for the shear viscosity and employing the previous ansatz for the relaxation time, 
we get 



G{T) 

c a 
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where G(T) and Q(T) are two dimensionless functions of the temperature rispec- 
tively for the gluon and the quark contribution. Due to the fact that the coupling 
constant g(T) in this model is a logarithmic decreasing function with the tempera- 
ture we have that the above functions approach their asymptotic limits 
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G(T -> oo) = rf s 4r(4)C(4)/(30vr 2 ) « 1.404 
Q(T -> oo) = d g (7/2) r(4)C(4)/(307r 2 ) « 2.763 



Therefore the analytic behavior oirj/s in Eq. (17) shows that the quasi-particle 



approach in the relaxation time approximation is able to reproduce asymptotically 



the correct pQCD limit. At very high temperatures, comparing Eqs. (16) and 
(17), we obtain the following values for the coefficients of the relaxation time C g = 



3.573 • 10-2, C q = (4/9)C s = 1.588 • 1(T 2 and a 2 = (/i*/T) 2 . In the T ^ oo limit, 



with these coefficients we obtain the same perturbative limit of Eq. ( 16 ) for both the 
hotQCD and Wuppertal-Budapest fits, even if the quickness with which this limit 
is reached can be quite different. Notice that in the asymptotic limit g{T)T — > oo, 
the same pQCD limit is recovered independently of the a 2 values, because in such a 
limit In^Ja^j gT) — > —IngT. 

Choosing //* = tcT as in Ref. |44j, the shear viscosity shows a minimum above 
T c for both fits and the value of r]/s remains close to the perturbative limit in fact 
we get a minimum (rj/s) m i n ~ 0.68 for the hotQCD fit and (r]/s) m i n ~ 0.66 for the 
Wuppertal-Budapest fit. This result however is in strong disagreement with recent 
lattice QCD calculations that show that the ratio 77/s ~ l/(47r) at temperatures 
close to the critical temperature, as well as with all the estimates of rj/s derived 
from heavy-ion collisions at RHIC energies within hydrodynamics [511 E2] and/or 
the parton transport model [53l [54] . 

Therefore, in a way similar to the g 2 parametrization, the value of the minimum 
in rj/s depends on the value of the a 2 coefficient in the logarithm. If we choose a 2 = 
25 (close to the value used for the other ansatz for r) while keeping the same values 
of the coefficients C g and C q in order to have the correct pQCD limit, we obtain that 
the minimum of 77/s ~ 1 / (Air) , in agreement with the current estimates in the 1 — 2 T c 
range from heavy-ion collisions phenomenology and first lattice QCD analyses. We 
see that the g A ansatz produces a significantly stronger and more rapid rise of 77/s 
(dashed lines in Fig|6]). There is however still a noticeable difference between 77/s 
in the WB case, that is always smaller than that calculated for the hotQCD fit. 
This is in agreement with the general slower convergence of the thermodynamics 
towards the perturbative gas limit in the WB case, that is reflected by the different 
behavior of the coupling constants, see left panel of Figj2j As shown in Figj6} we 
obtain that the position of the minimum for the dashed curves is very close to the 



one obtained in the previous parametrization for the relaxation time Eq.(14), and 



for increasing temperatures we have that rj/s is a monotonically increasing function 



for both fits. It asymptotically approaches the pQCD limit of Eq.(16), while for 



the parametrization of Eq.(14) we have that r]/s never reaches this limit. The 
hotQCD case gives a 4irr]/s ~ 10 (relatively close to the perturbative limit) already 
at T ~ 4T C hence predicting a nearly perturbative plasma in the initial stage created 
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Figure 7: Bulk viscosity to entropy density ratio (/s as a function of T/T c . The 
solid lines are for the r of Eq.(14), the dashed lines are for the r of Eq.(15). The red 
curves correspond to the hotQCD data, the black curves to the Wuppertal-Budapest 



ones. The dot-dashed curve corresponds to the NJL model with rj/s — (47r)~ 
the text for details. 



sec 



at the LHC The WB case instead shifts this limit at much higher temperatures and 
at T ~ 4T C we still have a 4irri/s ~ 4. The latter could be the reason behind a 
similar amount of observed elliptic flow at RHIC and LHC energies |37] 

To evaluate the bulk viscosity, Eq.(12), we have used the crude approximation 
of small relaxation time as for the shear viscosity. However, the bulk viscosity has 
a different origin and the processes involved could lead to a different r as discussed 
in Ref.[lH]. 

In Fig. [7] we show the bulk viscosity to entropy density ratio (/s as a function 



of T/Tc- As we can see from Eq. (12), the bulk viscosity is more sensitive to the 
equation of state than the shear viscosity and the term [Ei — T^f\ goes to p 2 / E for 
T ^> Tc; therefore, at high temperatures the behavior of the bulk viscosity depends 
on how C 2 — » 1/3. Again the difference can be traced back to the different quasi- 
particle masses, but in general the differences among the cases considered are much 
smaller than in the case of the shear viscosity. In general, there is a compensation 
between the behaviors of g(T) and C S (T) with a tendency to have a smaller (/s for 
the hotQCD case, that can be traced to the faster increase of C s towards 1/3. 

For comparison, in Fig. [7] we show the ratio (/s for the NJL model with the 
same set of parameters used for the speed of sound, see Fig. |4| In this calculation 
the relaxation time is chosen in such a way to keep the ratio rj/s constant and 
fixed to the lower bound rj/s = (47r) -1 . As we can see, the bulk viscosity for the 
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NJL model is much smaller by orders of magnitude than the one calculated using 
the quasi-particle model and it has a rapid decrease for T > Tq. This is a direct 
consequence of the fact that a large mass (m ~ 300—400 MeV) in NJL is related only 
to the chiral phase transition which modifies the speed of sound only for T < 1.2Tc 
and approaches the relativistic limit of C| = 1/3 just above T c , as shown in Fig[4j 
Note that such a large difference does not depend on the particular relaxation time 
used, because (/s for T > 1.2Tq is already several order of magnitudes smaller with 
resepect to the quasi-particle model. 



4 Susceptibilities 

We conclude our analysis by showing the predictions on quark number susceptibili- 
ties that we obtain in the quasi-particle model for the two different fits to the lattice 
data. Quark number susceptibilities are a very useful tool to understand the nature 
of the degrees of freedom in the vicinity of the QCD phase transition. They have 
been studied both on the lattice [2U H6] and using phenomenological models [I9j [50] 
The ones that we address here are defined in the following way: 



uu 
On 



Td 2 \nZ 



V dfj* 



s _ Td 2 \nZ 



The results published so far by the hotQCD collaboration on light quark suscepti- 
bilities [21] seem to indicate that quarks behave as almost massless, non-interacting 
quasiparticles immediately above Tc- This can be seen in Fig. [HJ where we plot 
the predictions of the quasi-particle model for light and strange quark number sus- 
ceptibilities, in comparison to the available lattice data from the hotQCD and the 
Wuppertal-Budapest collaborations (light quark number susceptibilities are not yet 
available in the latter case). 

Earlier works have calculated the quark susceptibilities using a quasi-particle 
model. In particular in Ref. [19] the g(T) of the quasiparticle model has been 
directly fitted to df^T), then predicting c 4 (T). Such a study has shown that, in 
principle, it is possible to account for the quark susceptibilities in a quasi-particle 
approach. Our study is however meant to check deeply the consistency of the present 
quasi-particle model not fitting the susceptibilities but calculating them from the fit 
to the energy density (see Fig. [2]). We can see from Fig. [8] that, for the hotQCD 
fit (dashed line) we underestimate the data. Indeed, the lattice data already reach 
the ideal gas limit for temperatures slightly above Tq, thus leaving little space 
for a thermal quark mass. In fact, such observable is very sensitive to the quark 
masses that we use as inputs. This apparent inconsistency between the description 
of QCD thermodynamics and susceptibilities within the quasi-particle model has 
been observed before 
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Apart from the comparison to the hotQCD data, our results show that we should 
expect for the WB case a smoother behavior of c^(T) across T c ; in fact we can see 
in Fig. [8] (left) that the susceptibilities corresponding to the Wuppertal-Budapest fit 
(solid line) are significantly lower than the hotQCD ones even at large temperatures, 
due to a difference of ~ 100 MeV in the corresponding quark masses (see Fig. [2J. 
Notice that the hotQCD lattice data shown in Fig. [8] are obtained for heavier-than- 
physical quark masses and at finite lattice spacing: they might change once the 
continuum limit is taken and physical quark masses are used in the simulations. 

In Fig. [8] we also show the results taking into account a finite strange quark 
mass (mo s = 150 MeV). This in principle makes the model more realistic. We have 
performed the same procedure of fit of the energy density as described in Section 2 
but for the non vanishing m Qs . For the pressure, interaction measure, gluon and ligth 
quark masses we get essentially the same result as in the vanishing m 0s case. We 
can see this also for C2 U (T) in Fig. [8] (left), where the mo s = 150 MeV case is shown 
by dashed lines and is quite similar to the solid lines corresponding to mo s = 0. 
Nevertheless, the effect is significant for the strange quark susceptibility \% as one 
would expect, but this increases the difference with respect to the available lattice 
QCD results for both WB and hotQCD. This is just another confirmation that the 
source of the discrepancy is the large value of the quark thermal masses. Therefore, if 
the trend of the present lattice data should be confirmed, in order for a quasi-particle 
description to hold, the thermal mass of quarks should be much lighter than the one 
we obtain from our fit. This would imply that the perturbative relationship Eq. 
Q between the quark and gluon thermal masses should not hold anymore, or some 
other dynamical mechanism should be included in the quasiparticle description. One 
candidate could be the Polyakov loop, whose dynamics can probably account for the 
gluon dynamics close to the phase transition. 

Another possibility could be to use the quark and gluon masses that come from 
the Hard Thermal Loop formalism. In this case, quark masses would be smaller 
than the ones we obtain from our fit, while gluons would be heavier as given by the 
following relations: 

In Fig. [9] we compare the susceptibilities that we obtain with the HTL relations 
(dash-dotted lines) to the pQCD ones from Fig. [8] (solid lines). In both cases, 
a fit to the lattice QCD energy density has been performed. It is clear that the 
results obtained by using the HTL ratio between quark and gluon thermal masses 
are significantly closer to the 1QCD data available for the hotQCD case, in particular 
for the strange susceptibility 
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A WB data 




T (MeV) T (MeV) 

Figure 8: Light (left panel) and strange (right panel) quark number susceptibilities 
as functions of the temperature. In both panels, the green lines correspond to the 
hotQCD fit, the black lines to the Wuppertal-Budapest one. The solid curves have 
been obtained by taking equal masses for the light and strange quarks, the dashed 
curves correspond to a heavier strange quark mass, with a shift m 0s = 150 MeV. 
The hotQCD data are taken from Ref. [21], the WB ones from Ref. [15] . 



5 Conclusions 

We have given an interpretation of recent lattice data for QCD thermodynamics in 
terms of quark and gluon quasi-particles with a thermal mass which we obtain from 
a fit to the lattice results. A thorough comparison between the results obtained from 
the fit to the hotQCD and Wuppertal-Budapest data is performed. We find that the 
differences between the results of the two collaborations are reflected in the differ- 
ences between the quasi-particle characteristics, such as masses and bag constant. 
In particular, the fit to the Wuppertal-Budapest data implies larger masses but with 
a weaker temperature dependence, which is associated to a weaker evolution of the 
running coupling. This is reflected on the dynamical properties of the quasi-particle 
medium. The shear viscosity over entropy ratio corresponding to the Wuppertal- 
Budapest fit shows a slower increase with increasing temperature with respect to 
the one corresponding to the hotQCD fit, thus favoring a hydrodynamical approach 
also in the LHC temperature regime. However we find that, once the quasi-particle 
properties are fixed to reproduce the energy density and/or the interaction measure, 
the predictions for quark number susceptibilities underestimate the presently avail- 
able lattice results for both the WB and hotQCD cases. This seems to indicate that 
the thermal quark masses that we obtain from the fits to QCD thermodynamics 
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Figure 9: Light (left panel) and strange (right panel) quark number susceptibilities 
as functions of the temperature. In both panels, the red lines correspond to the 
hotQCD fit, the black lines to the Wuppertal-Budapest one. The solid curves are 
the ones presented in Fig. |8j The dot-dashed curves are obtained by using HTL 
masses for quarks and gluons. The hotQCD data are taken from Ref. [23], the WB 
ones from Ref. [46J. 



are too large. This last result casts some doubts on the capability of the present 
quasi-particle model to correctly describe the inner dynamical structure of lattice 
QCD and hints at the necessity to further augment the dynamical structure of a 
quasiparticle approach in order to have a safe and reliable tool to understand the 
physics behind the lattice QCD data. A HTL approach suggests a smaller quark 
to gluon mass ratio m q /m g , which reduces the discrepancy with 1QCD data for the 
susceptibilities while mantaining the same level of agreement for the pressure and 
energy density. We notice that transport equation self-consistency associated to the 
present quasi-particle model can be easily derived in a way similar to Refs. f29l [32] . 
This will allow to go beyond the cascade approach that takes into account only 
collisions among massless partons [S31 EH ES] , including a field dynamics that can 
correctly reproduce the thermodynamics of 1QCD data. 
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